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Abstract 

We consider the extension of the standard model with an arbitrary 
number of ^7(1) gauge fields coupled to baryon-minus-lepton number 
and/or hypercharge. Under the assumption that A'^ from the LEPl 
experiment is an unlucky fluctuation, we find moderate evidence for 
the presence of such fields in the precision electroweak data. A rela- 
tively large range of the Higgs boson mass is allowed. We discuss the 
phenomenology of the extra ^7(1) fields. 

1 Introduction 

The present day data from the high-energy colliders like LEP and the Teva- 
tron show that almost all data are described by the standard model at the 
loop level. Therefore the extensions of the standard model (SM) tend to be 
strongly constrained. Typically such extensions will spoil the agreement with 
the data through a variety of effects, one of the most important of which is 
the appearance of flavor changing neutral currents. Even the most popular 
extension, the minimal supersymmetric extension of the standard model, has 
to finely tune a number of parameters. This leaves only one type of extensions 
that are safe, namely extensions with singlet particles. Since the discovery 
that neutrinos are massive, it is clear that singlet fermions play a role in 
nature. There are arguments from cosmology, that singlet scalars could be 
important [H [2], [3l HI [5] . In a recent analysis of the Higgs-search data from 
LEP2 it has been pointed out, that such singlet scalars may already have 
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been seen as a smeared-out Higgs boson [6J. Given this situation it is nat- 
ural to ask whether also singlet vector bosons can be present. Therefore we 
decided to study the most general renormalizable extension of the standard 
model containing extra gauge bosons, but no extra fermions or scalars. Since 
the mass of these extra bosons is put in without a Higgs mechanism, these 
can only be abelian vector bosons. Demanding the absence of anomalies in 
the gauge currents, the extra vector bosons can only couple to linear combi- 
nations of hypercharge {Y) and the difference of baryon and lepton quantum 
numbers [B-L). However they can have an arbitrary mixing with the stan- 
dard model hypercharge field. Through this mixing, these fields introduce 
small changes in the couplings of the Z-boson to matter. These couplings 
are suppressed by a factor m|/m|,. If the masses of the extra bosons are 
large enough such effects are allowed in the data. Of course if the couplings 
are infinitesimal, such an extension is always possible. The real question is 
whether the effects of the new physics can improve the agreement with the 
data. Even though almost all data are described by the SM |7j, there is the 
somewhat disturbing situation, that the overall fit to the standard model 
is not very good. The reason is the large difference in sin^ 6'^'^* from the 
forward-backward asymmetry of the bottom quarks and the measure- 
ment from the SLAG SLD experiment. No realistic model has so far been 
able to explain this difference. It has led some authors [8], IHl [lOl [HI [12], to 
reanalyze the data leaving out A^. In this case one can get a good fit to the 
data, however with an unphysical Higgs mass of 50 GeV. What we will show 
in this paper, is that we can get a good fit to the data for a physically allowed 
Higgs mass, if we allow for one or more extra U{1) fields. We must however 
emphasize that this works only in the reduced data set without A^. Using 
the full data set it appears that one has, close to the 95% GL, an unresolvable 
problem. Within the reduced data set we were able to find confidence level 
contour plots for the couplings of the extra vector bosons, as a function of 
the Higgs mass. One can actually allow for a larger range of the Higgs boson 
mass than in the SM. 

The outline of the paper is as follows. In section [21 we describe the 
models. In section [S] we discuss the Z-boson couplings in relation to the 
precision experiments. In section H] we discuss the phenomenology of a single 
Z'-boson. In section [5] we consider some interesting possibilities with more 
than one Z'-boson. In section [6] we recapitulate our conclusions. 
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2 The model 



There is a large class of models containing extra neutral vector bosons [13] . 
We limit ourselves to the simplest, so-called non-exotic extensions p^fT5l lT6l 
[T7] . In this type of extension the extra gauge bosons couple universally to the 
fermion generations and no exotic fermions are present. The model consists 
of the SM plus an extra number n of U{1) gauge fields. We call the U{1) fields 
C^, where i = 0,1, ... ,n. These fields can each have arbitrary couplings to 
the Y and to B-L. Other couplings are not possible, because these violate 
the renormalizability of the model through the anomaly. Besides this, we are 
concerned with the ordinary SU{2)l gauge fields W^. The pure gauge field 
part of the Lagrangian consists of a kinetic and a mass term. 
The most general form is: 

3 1 " 1 " 

a=l i,j=0 i,j=0 

where C^^, and are the field strength tensors for the U{l)s and SU{2)l 
fields respectively, and the % and mf^ are the coupling parameters. 

Through a linear transformation of fields {C C), we can bring the 
Lagrangian in the canonical form in the diagonal mass basis ( Zij —>■ 6ij, 
ihfj —>■ dijUif ). In order that the photon stays massless, one of the fields 
must be massless; we will call this field B^ = C^. The Lagrangian becomes 
therefore: 

1 ^ 1 1 " . . 1 " 

-^gauge = F^iu^fiu ~ '^B^yB^y — — C^j^C^j^ — — C^C^. (2) 

a=l 1=1 i=l 

Having the fields in the mass basis, the original couplings of the fields have 
changed. We have only one extra condition. This is that the B^ field must 
not couple to B-L, because otherwise the photon would have a B-L quantum 
number {Qb-l)- The fields have couplings gy to Y and g'^Q_j^ to B-L. 

The coupling to the fermions is then described through the minimal cou- 
pling 

(n n \ 

+ i^W^rl + iUsB, + Y,9yC;)y + 1^.93- lC^Qb-l ^. 
1=1 i=i J 

(3) 

Here s and c are sine and cosine of the weak mixing angle. The gauge boson 
coupling to the Higgs field originates from the same covariant derivative, but 
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the Higgs field has Qb-l = 0. After spontaneous symmetry breaking, we 
therefore find a mass matrix of the form: 



1=1 i=l 



We subsequently assume that rrii ^ m and diagonalize the mass matrix in 
perturbation theory in the parameters 



1 _ 

1^1 ~ c^T^l 



(5) 



In the following we will ignore effects of 0{1/ ^f). 

We recover the SM relation between the photon field [A^] and the La- 
grangian fields 

A^^sWl + cB^, (6) 
while for the Z-boson we find 

n j 

Z, = -cWl + sB,-J2 (7) 

i=i 

with the Z-boson mass given by: 

m| = — (1-ay), (8) 

where we defined 

".-t^. (9) 



For the Z'-bosons we find 



Z^-Cl + ^-^i-cWl + sB,), (10) 

with masses 

m% = ml (l + ^) . (11) 

With the relations above and Q = /2 + Y/2 we find for the coupling 
iip'jfj^g^ip of the Z-boson to the fermions: 

(3) 

4 = ^^(1 - (l - - cW) + 4rf g^-LaBy + 75) , (12) 



where we defined 



aBY 



i^yyy|^_ (13) 



The isospin quantum number r'^'^^ is +1 for neutrino and up-quark, —1 
for electron and down-quark. For the couphngs to the fermions we find 

-or / 1 \ 

4 = -\^9y f 1 - 4|Q| + 75 + ^(1 - m\s^ + 75) ] + gB-LQB-L. (14) 

For the Z-pole variables, ay and a^y are sufficient to describe the data, 
but for off-shell quantities also the direct Z'-exchange plays a role. We there- 
fore introduce the quantity 



75 (15) 
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In the case of a single extra U{1) boson one has the relation 

a|y = aYttB-L- (16) 

When one has more than one extra vector boson this relation is in general 
not true anymore. 



3 Z-boson Couplings and Fits to the Preci- 
sion Data 

In this section we compare the predictions of the model described above 
with the precision measurements of LEP and SLD at the Z-boson peak. We 
restrict our analysis to the case of a single Z'-boson. Since m^/ ^ rriz, 
the Z-peak observables are only sensitive to the modified Z-boson couplings. 
Effects from direct Z'-boson exchange are suppressed by a factor r|/m|,. 
The Feynman rule for the vertex coupling the Z-boson to a fermion / is given 
by 



i7^ {Vf + An,) , (17) 



where the vector and axial couplings Vf and Af are 

Vf = -^{l-aY-4\Q\{s'-aY)+4aBYr^'^QB-L), (18) 
Af = -^(1-ay). (19) 
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In the equations above, t'^'^^ = ±1 is the isospin third component, Q the 
electric charge and Qb-l is —1 for leptons and 1/3 for quarks. 

In perturbation theory, even at the tree-level, in order to make predictions 
for physical observables, it is necessary to fix the numerical value of the 
Lagrangian parameters g, s, and m in terms of the input parameter set a, 
mz and Gp- In the model under study, the fitting equations relating the 
Lagrangian parameters with the experimental data of the input parameter 
set are 

1 

Gf = (20) 
mz = — (1-ay). 

By solving the fitting equations, one finds that all the bare Lagrangian pa- 
rameters depend on ay 

g"^ = gl + kgtty, 

m = + kmay, (21) 

2 2 7 

S = - kstty, 

where the constant introduced above have the numerical values 

gl = 0.4322, kg = 0.5915, 

= 80.94 GeV, km = 55.38 GeV, (22) 
si = 0.2122, h = ^ = 0.2904. 

The theoretical predictions for the physical observables have to be ex- 
pressed in terms of the quantities g^, and (c^ = 1 — si). It is important 
to observe that, applying the Eqs. (1^ to the Z-boson coupling in Eq. (fT^ . 
one has to operate the replacement 

s - c ay s^ 2_ ^ ay. (23) 

In the case in which Q = —1, the quantity 1— 4s^ is close to zero and therefore 
sensitive to radiative corrections. This has been well studied within the SM. 
The main effect comes from the 7 — Z mixing terms from fermion loops and 
leads to a shift in the effective value of to s^.. ^ 0.2314, as measured in 
the fermion coupling to the Z-boson, which is the numerical value to use in 
comparison between the SM and the Z' model. The remaining uncertainty 
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in the prediction is of the order of the radiative corrections combined with 
Z' effects. Given the experimental errors affecting the measurements of the 
physical observables, the remaining uncertainties on the theoretical side will 
not affect the comparison with the data. 

Numerically we get the following relations 

^4' = i?4j^(l + 0.92ay-l.lW) (24) 
= (7^^|gj^(l-0.10ay + 2.28asy) (25) 

r^lz' = r^|sM(i + 0-iK-i-io«BY) (26) 

sin^ I ^, = sin^ Q^H^ | - 1 . lOay - a^y (27) 
mw\z, = miy|gj^(l + 0.72ay) (28) 



3.1 Analysis of the Z-pole data 

In order to test whether there is evidence in the data for this model we 
used the following precision measurements: the six sin^ 6^^^ measurements 
from LEPl and SLD, Ri, and Fz also from LEPl. Furthermore the 
niw measurements from LEP2 and the Tevatron were used. G^-. and 
mz are used as the input parameters of the model. Variation within la of 
TJit, Aajj^'* and ag do not affect significantly the following analysis. To see 
whether a given model, (i.e. a prescribed value of mn, ay and a^y) fits the 
data, we calculated the predicted values of the above mentioned variables and 
calculated for these values. We varied the parameters over the physical 
range mn > 115 GeV, ay > and gby arbitrary. 

We proceed as follows. We assume that the world is described by a known 
value of niH, ay and asy within the physically allowed range. On the basis 
of these values a prediction is made for each of the n physical observables Xj. 
We form the quantity 

n /^^measured _ ^predicted ^^ 

x'-i:^ '- (29) 

i=l i 

According to statistical theory this quantity should follow a Pn cumulative 
probability distribution 

roo yn/2-l -y/2 

Pn[x') = I iy\i^r^ (30) 

We can calculate the probability Pn(x^)- This number gives the prob- 
ability, that a random fluctuation in the data measurements would have a 
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probability smaller than the actually measured data. If this number is small, 
it means that the data are not well described by the assumed values of ttih, 
ay and a by- The 95% confidence level area in the space of rriH, ay and a By 
consists of the points with Pnix^) > 5%. It means that the points outside 
this range are ruled out with a "confidence of 95%" . 

In order to see whether any value of the parameters lies within 95% CL 
area, we looked for the lowest value of x^- The results are shown in Tab. 
([T]). For all models the lowest value of "was found at mn = 115 GeV, 
giving x^/d.o.f. ^ 19/11. This high value indicates a bad fit to the data, as 
was noticed before in [8], O [THl [HI US]. The presence of the extra U{1) field 
hardly improves the x^- This is in agreement with the result in [16] for the 
ttBy = case. Both the SM and the extended models are barely compatible 
with the data at the 95% CL. 

Table 1: Test of for 11 experimental tuh = 115 GeV. 



Model 




niH ay a by 


SM 
ay 


19.1 5.9% 
19.0 6.1% 
19.0 6.1% 


115 

115 3.3-10-5 

115 4.4-10-5 



There are two possibilities that a bad fit can arise. The first is that the 
model is simply not correct. The other possibility is that the model is correct, 
but that the data happen to contain an unlucky data point that spoils the 
fit. This last possibility can be tested for by removing one measurement from 
the data and seeing whether there is a large change in x^- If we find a very 
large change in by removing one point, we consider this as an indication 
that this point should not be used in the fit. 

We redid the analysis removing one of the data points in turn and the 
results are shown in Tab. ([2]) for the more significant changes. 

Only in the case of removing A'^ within the context of the extended 
models are we able to find a model that cannot be excluded within 68% CL. 
Confidence level contours are given in the Figs. ([IH2])- 

Altogether we consider this analysis as a mild statistical indication, that 
the point is an unlucky fluctuation and that the electroweak data are very 
well described by the SM with a light Higgs boson and one or more additional 
Z'-bosons. We notice that the largest part of the improvement comes from 
the presence of the ay term, the additional parameter a^y having a smaller 
effect. This allows for models of Z'-bosons all coupled to hypercharge only, 
which is consistent with renormalization. 



8 



Table 2: Test of for different sets of 10 data points. 



Excluded 


Model 


x' 


Ao(x') 


rriH 




ay 




(IBY 


^SLD 


SM 


14.8 


13.9 % 


134 


GeV 




— 




— 


^SLD 


ay 


14.8 


13.9 % 


134 


GeV 











^SLD 




14.1 


16.8 % 


115 


GeV 









-1.7-10-^ 




SM 


16.2 


9.5 % 


115 


GeV 










'-'had 


ay 


16.1 


9.6 % 


115 


GeV 


3.6 


■ 10- 


-5 




"^had 




14.0 


17.1 % 


194 


GeV 


8.9 


■ 10- 


-4 


-7.4 ■ 10"^ 


^fb 


SM 


13.1 


21.9 % 


115 


GeV 












ay 


9.8 


45.7 % 


115 


GeV 


3.0 


■ 10" 


-4 






0-Y, OjBY 


9.5 


48.8 % 


115 


GeV 


1.8 


• 10" 


-4 


1.6 ■ 10"^ 



3.2 Comparison with ofF-shell data 

Besides the Z-pole observables and the VT-boson mass, there are a number of 
low energy experiments that are in principle sensitive to the presence of a Z'- 
boson. Examples are atomic parity violation, the NuTeV experiment, M0ller 
scattering, and the {g — 2)^ experiment. We checked the sensitivity of these 
experiments to the Z' model. None of these experiments has a sensitivity to 
the ay and a^y variables, that is comparable to the sensitivity of the Z-pole 
observables. 

Also the LEP2 experiment is sensitive to the presence of a Z'-boson. A 
recent analysis of the LEP2 data for lepton production has claimed evidence 
for the presence of Z'-boson effects [T8l [T9] . Limits were put on the axial 
vector and vector couplings of the Z'-boson. Translated into our notation, 
the following conditions were found: 

9ay + 24:aBY + l6aB-L = (65 ± 30) ■ 10-^ (31) 
ay = (118± 158) • 10"^ (32) 

We see that the errors in these measurements are very large, so one cannot 
say that one is testing the model at the same level of precision as with 
the Z-pole data. The very large value of the couplings implied in the first 
formula is not very natural compared to the values we derived before. In 
principle, since we have the new parameter as-i, one could have an extra 
contribution coming from a second Z'-boson, coupled relatively strongly to 
B-L, but with little mixing. An alternative explanation could be that one is 
actually producing directly a spread-out Z'-boson. Such a model is presented 
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«5f(10-4) 




Figure 1: Significance plot for tlie parameter space (ay,aBy) for a fixed 
niH = 115 GeV. Tlie tliree contours contain a cumulative probability -Pio(x^) 
higher than 32%, 10% and 5% corresponding to the 68%, 90% and 95% CL. 
The point at the origin describes the SM, while the other two points represent 
the minimum of the distribution for the pure ay model and for the ay 
plus a BY model. 

in section 5. At the moment one should wait for an independent confirmation 
of the results in Ref. [iHl [I9]. 
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a) b) 




Figure 2: Contours of cumulative probability for Pio(x^) distribution for the 
parameter space {ay,aBy} at different Higgs masses, a) The different conn- 
tours describe 68% CL and from left to right masses of 115 GeV, 165 GeV 
and 215 GeV. The open line moves along the set of minima for the up to 
rriH = 217 GeV. b) The countours describe 95% CL and from left to right 
masses of 115 GeV, 300 GeV and 495 GeV. The open hne moves along the 
set of minima for the up to uih — 500 GeV 

4 Phenomenology of 

4.1 Renormalization group analysis 

In order to study the phenomenology of the Z'-boson one needs its mass, pro- 
duction cross sections, and decay branching ratios. These can be computed 
when coupling constants and mass are known. From the precision data we 
can put limits on the relative values of the B-L and Y couplings. But since 
the limits from the precision measurements are in the generic form 5f^/m|, 
we have no direct information on the mass itself. For low masses one can find 
limits from direct searches. At first sight there appears to be no way to give 
an upper limit to the mass of the Z'-boson. A very heavy boson might have 
a large coupling giving similar effects to a light boson with small couplings. 
Therefore, to get an upper limit, we need more theoretical input. Indeed 
it is the case that one cannot arbitrarily enlarge the couplings. Since the 
vector bosons are abelian, their effective couplings grow via renormalization 
group and will become infinite at the Landau pole. If one demands that the 
Landau pole is much larger than the electroweak scale, one can get an upper 
limit to the mass of the Z'-boson. We will give somewhat qualitative limits 
using the one-loop running of the couplings constants due to the fermions 
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and demand that the Landau pole hes beyond the Planck mass scale. In this 
case the new U{1) fields would not be too different from the SM gauge fields. 
The analysis is slightly complicated, due to the fact that the renormalization 
group mixes the different fields [Sni EH [22]. We first define the following 
auxiliary variables 



91 = 9-, 92 = g—- (33) 



The renormalization group equations become 



dt 
dgB-L 
dt 



9B-L {Sgl.L + ^gl + ^g2gB-L) , (34) 



^ = ^1 (%B-L + 10^2) + ^2(8^|_i + + 8g2gB-L) ■ 

Their solution is given algebraically by using the derived equations 

5 5-2 + 4:gB-L . . 
K — = constant, 

gB-Lgt 



5gl + 5(7? - 8g%_L 



glWi 
1 



constant, (35) 



— + lot = constant. 



By demanding that there is no Landau pole before the Planck mass we find 
the allowed region in the (72(0) versus 5'b-l(0) plane given in Fig. (JSj). 

Unfortunately, the constraints we get from this condition are not very 
strong. For example, if we take the a^y = case, one finds the limit 
(72(0) < 0.42. This corresponds to a limit 



2 

m|, < 0.32^. (36) 
ay 

So, even for the relatively large value of ay = 4 ■ 10~^, one finds a limit 
of mz' < 2.6 TeV. For smaller values of ay the limit gets correspondingly 
higher. Therefore, even though it is plausible that the Z'-boson is in the low 
TeV region and thereby within reach of the LHC, this is not quite guaranteed. 
Stronger upper limits on the Z'-boson mass can only be derived, when more 
theoretical information is put in, for instance by demanding some unification 
of coupling constants. Assuming Eq. (1311) to be true we can get better limits. 
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^2(0) 




gB-dO) 



Figure 3: Region of allowed coupling strength at the weak scale. 

due to the effect that in this case the ub-l term is the dominant one. Using 
similar reasoning as above one finds 

mz' < SOOtllo GeV. (37) 

In this case the Z'-boson could well be within reach of the Tevatron at high 
luminosity. 

4.2 Collider limits 

Nowadays the best limits on the presence of a Z'-boson come from the Teva- 
tron collider. The principle guiding the search is straighforward [23] . One 
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uses the decay of the Z'-boson into an electron-positron pair and looks for 
a peak in the invariant mass. In addition one can use information on the 
forward-backward asymmetry of the leptons. One therefore needs a predic- 
tion for the production cross section and for the branching ratio into leptons. 
Moreover, when one is working within a specific model, one can use the distri- 
bution of the leptons in the center of mass angle cos^*, leading to somewhat 
different limits for different models. This search has recently been made at 
the Tevatron [23] for the class of models discussed in Ref. [IT], but not for 
the specific models presented here. 

In the following we use the narrow width approximation and ignore in- 
terference with Z-boson and photon exchange. The production cross section 
times branching ratio (BR) is well described by the formul30: 

a{pp ^ Z' ^ e+e") = -^BR{Z' e+e") x 

{cuWui-s, Ml) + CdWd{s, Ml)) , (38) 



where 
and 



+ K\ (39) 



c, = + (40) 



In the above formulas the V and A' are the vector and axial-vector couplings 
of the quarks to the Z' as given in Eq. (IHj). Furthermore Wu and Wd are the 
luminosities for up-type quarks and down-type quarks, s the total energy 
squared of the collision in the center of mass system. 

The predicted branching ratio to electron-positron is given in our model 
approximately by 

^''-^-^- ^,nz'^fir '''' 

where we have neglected the effect of other than fermionic-pair decays. For 
generic massive fermionic-pair decay of a single Z' we have 







127rV 




m| 



2 

(v;' + A'') + 3^(r;'-A'') 



(42) 



*c.f. with Eqs. (3.8-9) from [13 
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At the parton level one can give simple formulas for the forward-backward 
asymmetry in the center of mass system. 



where the function Rf = V^j has the following values 

r -3-4^, for/ = e, 

-l-ll^' for/ = «, (44) 

For collisions the measured asymmetry is weighted with the quark lu- 
minosities. In collisions the overall asymmetry disappears and one must 
resort to rapidity dependent asymmetries. 



Now we address the question, whether the limits derived in section [3] can 
be used to make general constraints on the expectations at the Tevatron. 
Actually the range of parameters as implied in Fig. ^ is too large to give 
much of a prediction even at a 68% CL. Therefore to get an impression of 
a possible explicit phenomenology we present results with additional con- 
straints. Actually the analysis of the precision data has shown, that the 
biggest improvement comes from the introduction of the ay parameter, the 
additional a^y parameter giving a smaller improvement to the fit to the data. 
It is therefore reasonable to consider the phenomenology, assuming qb-l to 
be absent. This is a consistent condition under renormalization and forms 
an interesting class of models by itself. 

Since the properties of the Z'-boson are strongly correlated with the mass 
of the Higgs boson, we are interested in knowing what the most likely expec- 
tation tells us for different Higgs masses. Therefore we also consider the set 
of models which move along the minimization-line with increasing mu- 

In principle the model dependence is largely contained in the parameters 
c„ and Q, defined above. As argued in Ref. |T7] it would be useful to have 
lower limits on the Higgs-mass presented in the c„ — q plane. Unfortunately 
such a comparison has not been presented in the literature. Instead we 
proceeded estimating the lower limits for the allowed mass of the Z'-boson as 
a function of either ay or mH-, depending on the restrictions mentioned above. 
We calculated these ranges as follows: we used a LO program to predict the 
total cross section as a function of the coupling constant and the mass. Then, 
this cross section was normalized to the sequential Z'-boson, which has a 
lower limit of 850 GeV [23]. We then connected the lower bounds of m^/ 
within the two studied models. The derived bounds are of course somewhat 
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qualitative, because the sensitivity to the cos 6** distribution is not exactly 
modeled this way. A precise analysis would require taking into account a 
bidimensional distribution, including the angular one as well, whereby one 
cannot ignore detector effects [23]. However this needs a detailed simulation 
of the detector and comparison with the actual data, which is beyond the 
scope of this paper. To derive the upper limits we used the results from the 
previous section. 

In the analysis of the one-parameter model (a^y = 0), we limit ourselves 
to the 68% confidence interval 

0.8 • 10"^ < ay < 6.5 • 10-^ (45) 

which we used to scan the possible lower and upper limits for the Z'-boson 
as shown in Fig. (jl]). 

As an alternative we considered the best-fitted models for different values 
of rriH. In these models the parametrized equations used were 

ayiniH) = ("l^^ - 7 log (^) + 7.3 log^ (^)) ■ 10-^ (46) 

aBYimH) = (-42.2 + 411og(^)-6.11og^(^))■10-^ (47) 
being accurate in the interval 

115 GeV <mH< 500 GeV. (48) 

Here in order to extract the maxima of the (72(0) coupling domain, we pro- 
jected the best estimates of ay and a^y onto the ellipse of Fig. ([3]), obtaining 
thus a non-linear mapping tuh i — > g-y, which finally can be transformed, 
similarly as in Eq. ( 136|) . to a mz' upper limit. 

The result for this analysis is shown on Fig. . In order to avoid a wrong 
interpretation, we point out that the limits shown in this picture do not 
correspond to equal values of confidence level at each mn- They just provide 
the range for the best fit. Actually the confidence limit from the precision 
data fit ranges between 51% at ttih = 115 GeV and 95% at mn = 500 GeV. 

Even though the approach is somewhat limited, it suffices to retrieve 
three main features. First, the Tevatron limits a large part of the a priori 
allowed masses. Second, there is a range that could still lead to discovery 
of a Z'-boson at the Tevatron with more integrated luminosity. Third, we 
have not enough restrictions on the upper limit to guarantee discovery at the 
LHC. The reason is that ay can be quite small. However if the Higgs mass is 
larger than allowed within the standard model fits, then ay becomes larger 
and one should be able to find a Z'-boson at the LHC. For the Z' search at 
the LHC, the models we discussed are similar to other models and present 
no particular difficulties, therefore we refer to Refs. [Ml [25] . 
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Figure 4: Lower and upper bounds on the Z'-mass for models with gs-L — 0. 



mz'(TeV) 


3 




2 




1 




0.5 


150 200 250 300 350 400 450 500 



Figure 5: Lower and upper bounds on the Z'-mass for the best model values 
{ay, Ofiy) as a function of uih- 
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5 Multiple Z^'-bosons 



The formalism as presented in section [2] leaves open the possibility of having 
more than one extra vector boson. From the precision data there is no 
possibility to get information about the actual number of Z'-bosons, because 
one measures only the parameters ay and a^y, that are formed by sums 
over all existing vector bosons. If one adds only a finite number of vector 
bosons there will be more Z'-bosons with different sets of couplings and 
the renormalization group running becomes more complicated. However no 
essentially new aspects appear. The situation becomes more interesting when 
one allows for an infinite number of extra fields. In particular one can also 
allow for vector fields moving in more than four dimensions. However for 
general couplings this will lead to a non-renormalizable theory, since higher 
dimensional gauge fields coupled to fermions form operators of a dimension 
higher than four. Here the comparison with the data from section [3] gives a 
useful hint. It was found that the major improvement in comparison with the 
data comes from the introduction of the ay term. The additional qby term 
gives a smaller improvement in x^. It is therefore natural to consider a model 
where all U{1) fields couple to hypercharge only. In the mass-basis, that we 
used so far, we then have a large number of fields all coupled to hypercharge 
with a certain strength. In this case however it maybe advantageous to 
study the model in the hypercharge basis. One then has one field coupled to 
hypercharge and further fields coupled to nothing. These extra fields make 
their presence known only through the mass mixing with the hypercharge 
field. Since a mass term for a vector boson has a mass dimension ci — 2 as an 
operator, one can allow for such vector fields as long as ci < 6. As an example 
we therefore take an abelian sector with a 4-dimensional field i?^, coupled to 
hypercharge and a (i-dimensional field coupled to nothing. We allow for 
a 4-dimensional mass term, a mixing mass term and a (i-dimensional mass 
term. 

The Lagrangian becomes 



i=l 

We consider the case that the fields move in a fiat open space. By first 
compactifying the higher dimensions and subsequently taking the continuum 



•gauge 




(49) 
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limit one can derive a hyper charge-boson propagator with a nontrivial Kallen- 
Lehmann spectral density. The analysis follows the treatment of the Higgs- 
singlet mixing in Ref. [6]. For a detailed derivation we refer to this paper. 
Here we only give the result. 

The hypercharge propagator becomes of the form 



d-6 



D^.'^iq') = V + ^^l,\q' + m')^ \ . (50) 

The masses M, m and /iihd are the free parameters of the model. The scale 
/iihd stands for low-to-high-dimensional mixing mass and measures the mixing 
of the high(=(i) dimensional vector and the low(=4) dimensional vector. The 
propagator contains a particle peak and a continuum. In order to guarantee 
that the particle peak is at mass zero, which is necessary to have a massless 
photon in the theory, we must take: 

M^m^""^ = (51) 

We now consider the simple cases of integer dimensions 4, 5 and 6. In the 
four dimensional case we should recover the original model. This is indeed 
what happens. One finds a single Z'-boson, with 

m|, = + m^ (52) 

and ^ 

ay = sin^ Ow ^ , ..^ , ^. (53) 

Next we consider the case ci = 5. The propagator becomes of the form: 



-1 



(r) = 5/.. [q +M'- mM\q' + m')-2j . (54) 
This corresponds to a Kallen-Lehmann spectral density: 

2rn? 6{s — rn?) mM"^ {s — rn?)^ 

P^^> = 2m^ + M^ ^'^^ ^ n {s - m^) {s - M'^)'' + ' ^^^^ 

We therefore find one massless particle and a massive continuum. After 
spontaneous symmetry breaking the massless excitation becomes the photon. 
The massive fields form a spread-out hypercharge-coupled Z'-boson. As a 
consequence we find: 

2, ds (2m2 + M2)M2(s-m2)l ^ ^ 

ay = sm Ow / 7^ 7 o^T^ jirn- 56 
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Finally we consider the case d = 6. This case is special, since it corre- 
sponds to the limiting dimension, where the theory is still renormalizable. 
Using a limiting procedure around d — 6 the propagator can be written as: 



I 2 I 2 N "1 ~1 

2 , ^2 2 1 q +m'' 



(57) 



The spectrum has a massless pole when 



2 

m 



+ /.f,, log ^ = 0. (58) 
The corresponding Kallen-Lehmann spectral density is: 

p{s) = J^^ , 5{s) + e{s - m") ^^ii^-^ . (59) 

+ \s- /.L log(^)] + TT^/^L 

Correspondingly one has: 

ay = sin^ dw r ds"^^-^^ . (60) 

[^-<dlog(^)] +7rVhd 

Since such a spread-out Z'-boson has no direct mass peak, its detection 
could be quite difficult at hadron colliders. Only an analysis together with 
an experimental detector simulation can give a reliable answer about the 
limitations in detection here. An important role could be played by the 
correlations for the outgoing leptons. At the Tevatron such a correlation can 
be measured directly, but at the LHC a rapidity dependent analysis will be 
needed. It seems even possible, that such a Z' could have been produced at 
LEP2, but has been overlooked because of its rather non-specific signature. 
A re-analysis of the data at LEP2 with this possibility in mind might be 
useful. 



6 Conclusion 

We studied the class of models of multiple U{1) fields coupled only to linear 
combinations of hypercharge and baryon-minus-lepton number, the so-called 
non-exotic U{1) fields. We took a careful look at the Z-pole precision data 
and argued the case, agreeing with some considerations in the literature, that 
one could reasonably analyze the data without the point. In the sub- 
sequent analysis wc found a moderate indication for the existence of extra 
U{1) fields. We studied the phenomenology of the extra Z'-bosons, conclud- 
ing that it is possible that the Tevatron could find the Z'-boson and likely. 
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but not entirely certain, that the LHC will find it. Furthermore we pointed 
out the interesting, but somewhat disturbing possibility, that the U{1) fields 
come from higher dimensions, leading to a spread-out Z'-boson, whose study 
might be severely limited at a hadron collider. If the energy would be high 
enough to produce it, a electron-positron collider could study such a spread- 
out signal without severe difficulties. 
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